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Results on charged particle production in p+p, d+Au and Au+Au collisions at RHIC energies 
(\/sjvjv = 19.6 to 200 GeV) are presented. The data exhibit remarkable, and simple, scaling 
behaviors, the most prominent of which are discussed. 



1 Introduction 

In this paper, we report on the bulk features of particle production, focusing on results from 
the PHOBOS experiment at the Relativistic Heavy- Ion Collider (RHIC). A primary strength of 
PHOBOS is the ability to measure the total charged multiplicity, with practically no extrapo- 
lations necessary, due to the ~ 47r acceptance of the detector. The specific detectors that allow 
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Figure 1: Excitation function of midrapidity (|r;| < 
1) charged particle yields per participant pair, 
dNch/dri/{\{Nj,art)), from PHOBOS at RHIC en- 
ergies for the 6% most central Au+Au collisions 
(solid circles). Also shown is heavy-ion data from 
lower energies and, for a baseline comparison, the 
midrapidity (jryl < 1) yield of charged particles from 
inelastic p(p) -I- p collisions (open diamonds). 
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Figure 2: Centrality, or (Npart), dependence of the 
charged particle multiplicity per participant pair at 
midrapidity for Au-|-Au collisions at 200 and 19.6 
GeV divided by the corresponding value obtained 
in inelastic p(p) + p collisions. The inset has the 
same axis labels as the main figure, with curved 
lines representing expectation from pure binary col- 
lision (Ncou) scaling and the horizontal dashed line 
corresponding to pure participant (Npart) scaling. 



the measurement of charged particle multipUcity {dNch/dr]) over the broad pseudorapidity range 
{\7]\ < 5.4) are the midrapidity Vertex detector (|?y| < 1), the centrally located Octagon barrel 
(|r/| < 3.2), and a series of six Ring counters (3.1 < |?7| < 5.4), all of which closely surround a 
thin Be beam-pipe. These detectors are silicon pad sensors, see Ref.'^ for details. 



2 Bulk Charged Particle Production 

A surprising discovery at RHIC is the existence of simple scaling behaviors observed in bulk 
charged particle production (i.e. integrated over p^ and particle species). Some of these features 
are similar to those seen at lower energies, as well as in simpler systems, including p(p) + p and 
e"'"-|-e~ collisions. These scaling properties are found to exist both as a function of the collision 
geometry as well as in the rest-frame of one of the colliding nuclei. In the following discussion, 
it is important to note that the measured bulk charged particle multiplicities are completely 
dominated by the emission of low p^ (< 1-5 GeV/c) particles. 

Determining the centrality (or impact parameter b) of a heavy-ion collision is extremely im- 
portant to provide a geometrical scale for use in the studies of underlying collision dynamics. The 
event centrality is characterized by charged particle multiplicities measured in various regions 
of phase-space. Comprehensive Monte Carlo (MC) simulations of these signals, that include 
Glauber model calculations of the collision geometry, allow the estimation of Npart, the number 
of participating nucleons in the collision 1=1^, for a selected class of events. The most central 
(b ^ 0) collisions will have the largest number of participants with the obvious upper limit of 
394 for a "perfectly central" Au-|-Au collision where all the nucleons interact. For a particular 
collision geometry, the MC simulation also allows calculation of the number of binary collisions, 
NcoU, which provides an expected baseline scaling for large- momentum transfer processes. 

One of the first measurements at RHIC was the multiplicity of charged particles, dNch/drj 
(pseudorapidity density) , produced in head-on collisions of Au-|- Au nuclei and emitted at midra- 
pidity, defined here as |?7| < 1. This and subsequent results are shown in Fig. ^ where the 
midrapidity dNch/drj is divided by the average number of participating nucleon pairs, {^{Npart)), 
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Figure 3: Centrality, or (Npart), dependence of the 
ratio of midrapidity charged particle multiphcity 
per participant pair for Au+Au coUisions measured 
at y/sNN = 200 GeV to that obtained at lower en- 
ergies y/SNN = 130 and 19.6 GeV. 
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Figure 4: Comparison of total charged particle mul- 
tiplicities per participant pair in p-(-p, d-|-Au and 
Au+Au collisions versus centrality at y^sjvjv = 200 
GcV. Also shown (dashed line) is the value from 
elementary e^-|-e~ collisions. 



in order to compare directly to results from p(p) -|- p collisions. What is most notable in the 
heavy-ion data of Fig. ^ is the apparent linear logarithmic growth of the midrapidity parti- 
cle production per participant pair with collision energy (-y/J/vJv)- This multiplicity in central 
Au-|-Au collisions is larger, however, than seen in inelastic p(p) -|- p collisions at similar energy, 
indicated by the open diamonds'*^ in Fig. ^ 

The centrality dependence of midrapidity multiplicity per participant pair in Au+Au col- 
lisions at ^SNN = 200 and 19.6 GeV is shown in Fig. 121 where the heavy-ion data has been 
normalized by the corresponding multiplicity in inelastic p(p) + p collisions. The most dramatic 
feature of the Au+Au centrality dependence is the similarity of the collision geometry scaling 
at both energies, a scaling which follows a strong A^pa^i-like dependence. Despite a factor of 10 
increase in collision energy from ^/sJ^ = 19.6 to 200 GeV, both data sets are consistent with 
only a small (~ 13%) fraction of binary collision (Ncou) scalingl^. 

Although there is a slight increase in midrapidity multiplicity per participant pair in Au+Au 
collisions as a function of centrality (see Fig. [21) , we find this increase exhibits a simple scaling at 
RHIC. This scaling is shown in Fig.|31 where the ratio of 200/130 and 200/19.6 GeV midrapidity 
dNch/dr]/{^ (Npart)) is given as function of centrality. Within errors, the ratio is independent of 
centrality for the most central 40% of the cross section and yields a simple scale factor R200/130 
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Expanding to the full phase space, we find the midrapidity increase in dNch/dri/{\ {Npart)) 
with centrality is remarkably compensated by a corresponding decrease in yield at high 77 such 
that the total charged particle production per participant pair is independent of centrality. This 
can be seen in Fig. |3] for both Au+Au collisions as well as for d+Au collisions. We also observe 
that the Au+Au and d+Au data give different values of Nch/{\{Npart)), but they are different 
in a very interesting way. For the same collision energy, the Au+Au data matches that seen in 
e++e~ collisions and the d+Au data matches that seen in p(p) + p collisions. 

The last scaling feature we examine is evident when particle production is viewed in the rest 
frame of one of the colliding nuclei. This is approximately achieved by effectively transforming 
measured dNch/dr] multiplicities to the rest frame of the beam with rapidity ybeam-, using dN/drj' 
where rj' = \rj\ — ybeam- In the late 1960's it was hypothesized that in hadronic collisions, 
at high enough collision energy, the charged particle yield (d^N/dy'dpy) reaches a limiting 
value and becomes independent of energy when viewed near beam rapidity^. This "limiting 
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fragmentation" hypothesis works in many elementary systems including p(p) + p as shown in 
the dNch/drj multiplicities in Fig. IHI We find the same features in the Au+Au heavy-ion data 
at RHICl^ (see Fig. |H1). The same "limiting fragmentation" features are also seen in dN^h/dr] 
data from d+Au collisions . Most recently, PHOBOS has additionally observed the same type 
of scaling in the magnitude of elliptic flow, V2{rj'), for heavy-ion collisions'^. 

3 Conclusions 

We report here a few of the bulk properties of charged particle production seen in relativistic 
heavy-ion collisions at RHIC energies. We find that much of the data can be described in terms 
of simple scaling behaviors. These observations are not yet understood, but could suggest the 
existence of strong global constraints or some form of universality in the bulk hadron production 
mechanism. These behaviors must be understood before a complete picture of the properties of 
the high energy-density matter created at RHIC can be definitively determined. 
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